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CN i ABSTRACT 

We present photometry on 23 Jupiter Family Comets (JFCs) observed at large helio- 
' C"| " centric distance, primarily using the 2.5m Isaac Newton Telescope (INT). Snap-shot 

O |. images were taken of 17 comets, of which 5 were not detected, 3 were active and 9 were 

I ■ unresolved and apparently inactive. These include 103P/Hartley 2, the target of the 

NASA Deep Impact extended mission, EPOXI. For 6 comets we obtained time-series 
photometry and use this to constrain the shape and rotation period of these nuclei. The 
^ ' data are not of sufficient quantity or quality to measure precise rotation periods, but 

the time-series do allow us to measure accurate effective radii and surface colours. Of 
the comets observed over an extended period, 40P/Valsala 1, 47P/ Ashbrook- Jackson 
and P/2004 H2 (Larsen) showed faint activity which limited the study of the nucleus. 
^j* Light-curves for 94P/Russell 4 and 121P/Shoemaker-Holt 2 reveal rotation periods 

of around 33 and 10 hours respectively, although in both cases these are not unique 
solutions. 94P was observed to have a large range in magnitudes implying that it is 
one of the most elongated nuclei known, with an axial ratio a/b > 3. 36P/ Whipple 
was observed at 5 different epochs, with the INT and ESO's 3.6m NTT, primarily 
in an attempt to confirm the preliminary short rotation period apparent in the first 
data set. The combined data set shows that the rotation period is actually longer 
than 24 hours. A measurement of the phase function of 36P's nucleus gives a rela- 
tively steep /3 = 0.060 ± 0.019. Finally, we discuss the distribution of surface colours 
observed in JFC nuclei, and show that it is possible to trace the evolution of colours 
from the Kuiper Belt Object (KBO) population to the JFC population by applying a 
'de-reddening' function to the KBO colour distribution. 
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1 INTRODUCTION 



The study of short period comets has recently benefited 
from the successful spac ecraft missions to 19P/Borrelly 
dSoderblom fc et aLlliooj). 8 IP /Wild 2 (Brownlee fc et all 



120041 ) and 9P/Tempel 1 (|A Hearn fc et alll2005t ). which pro- 
vided a wealth of information on these nuclei. The Rosetta 
space-craft is on its way to 67P/Churyumov-Gerasimenko 
(|Schwehm fc SchuhJ Il998f ) and will provide another great 



*Based on observations collected at the European Southern 
Observatory, Chile (proposals: ESO No. 74.C-0125, 278.C-5040), 
and at the Isaac Newton Group of telescopes on La Palma (pro- 
posals: I/2005A/11, I/2006A/7). 

fE-mail: csnodgra@eso.org 

JNow at: NASA Jet Propulsion Laboratory, Science Division, 
4800 Oak Grove Drive, Pasadena, CA 91101, USA. 



leap by studying the nucleus in detail while orbiting around 
it. However, these missions provide great detail on only indi- 
vidual nuclei, and have themselves demonstrated the great 
differences between the comets visited. We therefore wish to 
study a far larger population of cometary nuclei, in order to 
better understand what they have in common, and which 
features vary between them. 

Sending space missions to a larger fraction of the 
cometary population is clearly impractical, and too expen- 
sive to seriously consider. Instead, we seek to measure the 
bulk properties of many nuclei using ground based telescopic 
observations. Although the activity of comets masks the nu- 
cleus when they are bright and near the Sun, it is possible to 
observe Jupiter Family Comets (JFCs) when they are closer 
to aphelion (and apparently inactive) using 2-4m class tele- 
scopes. 
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A number of groups have used this method to take 
'snap-shot' observations of distant nuclei, allowing mea- 
surement of their brig htnesses and ther e fore, assuming 
an a l bedo, their s i zes (lLowrv et alj 120031 ; iLicandro et al.l 
l200d ; iMeech et ail I2004T ). We have expanded upon this 
by taking time-series observations of JFC nuclei; by ob- 
taining light-curves, constraints on the size, shape and 
bulk density of the nuclei can be measured. We also 
take multi-filter photometry in order to study the sur- 
face colou rs of nuclei. The first results of t his wor k are 
given by Snodgrass. Fit zsimmons fc Lowrvl |2005bl ) and 



ISnodgrass. Lowrv fc Fitzsimmons! l|20od ) (hereafter papers 
I & II). 

In this paper we present the results of observations (or 
attempted observations) of 23 comets. Section[3]gives the re- 
sults of further snap-shots obtained on comets that were not 
followed over an extended period, while section [4] gives the 
results of the time-series observations. We observed comet 
36P/ Whipple (hereafter 36P) on three separate occasions, 
and the results from this large data set are described in sec- 
tion Section HJ] discusses our findings in terms of results 
for other comets, draws conclusions about the population 
of JFCs as a whole, and compares JFCs with Kuiper Belt 
Objects (KBOs), the supposed parent population for JFCs. 
Finally, section [7] summarises all of our results. 



2 OBSERVATIONS AND DATA REDUCTION 

The data presented were taken using three different tele- 
scopes over 5 observing runs of varying lengths. A log of all 
the observations is given in table[2] The majority of the data 
were taken with the 2.5m Isaac Newton Telescope (INT) at 
the Roque de Los Muchachos observatory on the island of 
La Palma. 

The INT was used during two runs. The first was in 
July 2005, on seven nights (July 1st - 7th) around the time 
of the Deep Impact mission's encounter with 9P/Tempel 1. 
The telescope was being used to follow the evolution of the 
dust coma and tails of this comet before a nd after the im- 
pact (|Lowrv et al.ll2005l ; Ijones et all 120051 ) . but was avail- 
able during the later half of each night after 9P had set. 
Conditions were excellent and all seven nights were pho- 
tometric; time-series data were collected on three comets 
[40P/Vaisala 1, 94P/Russell 4 and P/2004 H2 (Larsen)] and 
snap-shot observations were made of 7 others. The telescope 
was used again between the 27th of February and the 2nd of 
March 2006, during which the first two nights were lost due 
to severe winter weather. The third night was clear but non- 
photometric, while the 4th night was photometric with peri- 
ods of excellent seeing (the FW HM of the image PSF varied 
between 0.8 and 2.7", with a median of 1.1"). Time series 
observations of comets 36P/Whipple and 121P/Shoemaker- 
Holt 2 were obtained, and snap-shots were taken of 12 other 
objects. 

The imager on the INT is the Wide Field Camera 
(WFC) which is mounted at the primary focus. The WFC 
is a mosaic of four thinned EEV 2048x4096 pixel CCDs, 
which has a pixel scale of 0.33" per pixel and a total field of 
view of 34' on a side. The wide field nature of the camera 
is not necessary for this work, and only one element of the 
multi-extension fits file is used to minimise the reduction 



workload. In the 2006 data this mostly corresponds to the 
central CCD 4, while the comets were centred on CCD 3 in 
the 2005 data as the telescope pointing had been offset to 
make use of the full width of the WFC mosaic in observing 
9P's tails. Even the single CCD has a field of view of 11.5 x 23 
arcmin 2 , which gives the advantage of having a constant set 
of comparison stars over two nights for all but the fastest 
moving comets, allowing light-curves to be fully differen- 
tial with no uncertainties from night to night calibration. 
Also, by a fortuitous coincidence, comets 44P/Reinmuth 2 
and 47P/ Ashbrook- Jackson were close together on the sky 
during the February /March 2006 run, and the wide field of 
the WFC allowed both to be observed simultaneously giving 
'two for the price of one'. 44P fell on CCD 1 and 47P fell on 
CCD 2. The filters used were the Harris V and Sloan r' and 
i! . The Sloan filters have the advantage of having sharp cut- 
on and -offs with little overlap in sensitivity at wavelengths 
between the bands, although further reduction steps are re- 
quired as the i'-band has significant fringing which had to 
be subtracted from the data manually for each frame, as 
pipelines or automated routines failed to measure the fringe 
brightness correctly. 

We used the 3.6m New Technology Telescope (NTT) at 
the European Southern Observatory's La Silla site on the 
nights of the 5th to 7th March 2005. This telescope, the 
EMMI instrument used, the run and the associated data 
reduction are described in paper II, where the majority of 
results from it were published. The data taken on 36P during 
that run were not included in paper II as it is more sensibly 
grouped with the other observations of this comet described 
here. The NTT and EMMI were also used in the 3rd time- 
series on 36P, taken on the nights of the 26th and 27th 
of February 2007. The conditions were clear, with seeing 
measured in the images between 0.6 and 1.3", with a median 
of 0.8" FWHM. 

The fifth set of observations was carried out using the 
2.0m robotic Faulkes Telescope North (FTN). We used ap- 
proximately 1 hour on the 31st of May 2006 to perform 
follow up observations on comet 121P/Shoemaker-Holt 2, 
which INT observations had shown to be brighter than ex- 
pected and showing evidence of having only recently ceased 
out-gassing (see section |4~4)| . to search for any sign of con- 
tinuing activity. The comet was imaged in the Bessell R- 
band using the DillCam CCD imager, which has a single 
2048x2048 pixel CCD with a pixel scale of 0.278 arcsec per 
pixel in the default 2x2 binning mode, giving a 4.6x4.6 
arcmin 2 field of view. Frames were also taken of the comet 
through the Bessell 1^-band, and of standard stars through 
both bands, to allow for photometric calibration. 

The reduction for the NTT, FTN and 2006 INT data 
and photometry on all frames were performed using stan- 
dard IRAF tasks as described in papers I & II. The 2005 
INT data were reduced (along with the 9P data taken dur- 
ing the run) using the Wide Field Survey (WFS) pipelin^; 
inspection of data from each run showed that both reduction 
methods gave equally good results. 

Differential photometry provides a precise measurement 
of how the brightness of the comet varies with time by com- 
paring it to the brightness of nearby field stars, which can 



1 http : //www. ast . cam. ac .uk/~wf csur/ 



Observations of 23 distant JFCs 3 



Table 1. Log of all observations. 



Comet 


UT Date 


Telescope 


r H [AU£ 


A [AU] 


a [deg.] 


7V cxp x Filter^ 




Appearance c 


_ ,/ 
8 pd 


01/07/05 


INT 


7.42 7 


6.76 


6.3 


13xr' 


95 


S 


36P 


05/03/05 


NTT 


4.08° 


3.39 


11.0 


5x_R, V, I 


200 


S* 




06/03/05 


NTT 


4.08° 


3.41 


11.2 


17xR, V, I 


200 


S* 




01/03/06 


INT 


4.78° 


3.79 


1.0 


63xr', lOxV, 10xi' 


150 


S 




02/03/06 


INT 


4.78° 


3.79 


1.0 


55xr', 5xV, 5xi' 


75 


S 




12/01/07 


NTT 


5.14° 


4.94 


11.0 


6x_R, 2xV 


220 


S 




26/02/07 


NTT 


5.17° 


4.35 


6.7 


84x_R, V 


120 


S 




27/02/07 


NTT 


5.17° 


4.34 


6.5 


8xR 


120 


S 




18/07/07 


NTT 


5.23° 


5.45 


10.7 


5xR, 2xV 


100 


S 


40P 


01/07/05 


INT 


4.57° 


3.68 


7.0 


8xr' 


85 


S* 




02/07/05 


INT 


4.57° 


3.68 


6.8 


15xr', V, i' 


85 


S* 




03/07/05 


INT 


4.58° 


3.68 


6.6 


7xr', V, i' 


85 


S* 




04/07/05 


INT 


4.58° 


3.68 


6.4 


20xr', V, i' 


85 


S* 




05/07/05 


INT 


4.59° 


3.67 


6.2 


17xr', V, i' 


85 


S* 




06/07/05 


INT 


4.59° 


3.67 


6.0 


20xr', 2xV, 2xi' 


85 


S* 




07/07/05 


INT 


4.60° 


3.67 


5.7 


18xr', 2xV, 2xi' 


85 


S* 


44P 


01/03/06 


INT 


4.51 7 


3.84 


10.1 


2xr' 


300 


F 




02/03/06 


INT 


4.50 7 


3.83 


10.0 


3xr', 3xV, 3xt' 


170 


S 


47P 


01/03/06 


INT 


5.11 7 


4.45 


8.9 


2xr' 


300 


S* 




02/03/06 


INT 


5.11 7 


4.43 


8.8 


3xr', 3xV, 3xi' 


170 


s* 


56P 


01/03/06 


INT 


3.83° 


3.14 


11.8 


3xr', 2xV 


520 


A 




02/03/06 


INT 


3.83° 


3.15 


12.0 


5xr', V, i' 


260 


A 


70P 


01/07/05 


INT 


4.84 7 


4.26 


10.6 


6xr' 


295 


F 


72P 


02/03/06 


INT 


3.28° 


2.35 


7.1 


5xr' 


50 


N 


75P 


02/07/05 


INT 


4.95 7 


3.94 


1.7 


3xr' 


75 


N 


78P 


02/03/06 


INT 


3.84° 


2.86 


1.8 


3xr' 


60 


A 


94P 


04/07/05 


INT 


4.14° 


3.19 


5.6 


14xr', V, i' 


75 


S 




05/07/05 


INT 


4.14° 


3.18 


5.3 


20xr', 3xV, 2xi' 


75 


S 




06/07/05 


INT 


4.14° 


3.18 


5.1 


26xr', 2xV, 2xi' 


75 


s 




07/07/05 


INT 


4.15° 


3.18 


4.9 


23xr', 2xV, 2xi' 


75 


s 


103P 


01/03/06 


INT 


5.03° 


4.30 


8.3 


2xr', V, i' 


200 


A 




02/03/06 


INT 


5.03° 


4.29 


8.1 


r', V, i' 


110 


A 


104P 


05/03/05 


NTT 


3.06° 


2.31 


13.8 


2xR 


80 


N 


114P 


04/07/05 


INT 


3.75 7 


2.95 


10.9 


8xr' 


90 


S 


120P 


02/03/06 


INT 


3.89° 


3.05 


8.8 


5xr' 


115 


F 


121P 


01/03/06 


INT 


3.92° 


3.43 


13.5 


4xr', V, V 


400 


S 




02/03/06 


INT 


3.93° 


3.42 


13.4 


18xr', V, i' 


215 


S 




31/05/06 


FTN 


4.18° 


3.34 


8.6 


UxR, 3xV 


60 


s 


131P 


02/03/06 


INT 


3.48° 


2.65 


10.1 


3xr', 3xV 


105 


s 


135P 


02/03/06 


INT 


3.55 7 


2.65 


7.5 


5xr' 


80 


N 


160P 


02/03/06 


INT 


3.98° 


3.41 


12.5 


3xr' 


140 


F 


P/1995 Al 


02/07/05 


INT 


5.52 7 


4.71 


6.9 


8xr' 


135 


S 


2006 BZS^ 


02/03/06 


INT 


2.32 7 


1.35 


7.1 


15xr', 3xV, 3xi' 


10 


S 


P/2004 H2 


01/07/05 


INT 


3.70° 


3.06 


13.6 


lAxr' 


210 


s* 




02/07/05 


INT 


3.70° 


3.06 


13.4 


19xr', V, i' 


210 


s* 




03/07/05 


INT 


3.71° 


3.05 


13.3 


14xr', V, i' 


210 


s* 




04/07/05 


INT 


3.71° 


3.04 


13.1 


16xr', 2xV, i' 


210 


s* 




05/07/05 


INT 


3.71° 


3.03 


12.9 


14xr', 4xV, 2xi' 


210 


s* 




06/07/05 


INT 


3.72° 


3.02 


12.7 


10xr', 6xV, 3xi' 


210 


s* 




07/07/05 


INT 


3.72° 


3.02 


12.6 


llxr', 4xV, 2xi' 


210 


s* 


P/2004 H3 


05/07/05 


INT 


3.71° 


2.97 


12.0 


6xr' 


90 


F 


P/2004 Tl 


02/03/06 


INT 


3.82° 


2.86 


3.9 


5xr' 


60 


N 


P/2001 X2 


02/07/05 


INT 


5.02 7 


4.02 


2.4 


2xr' 


80 


N 



a Superscripts / and O refer to whether the comet is inbound (pre-perihelion) or outbound (post-perihelion). The variations in R^, A 
and a over the course of any one night were smaller than the significance quoted here. 

^ The number of exposures -/V cxp in each filter in the final data set: frames rejected during reduction/analysis are not counted. 

c Appearance: S - Stellar, A - Active, N - Not detected, F — detected but too faint to determine whether or not activity was present: 

these are assumed to be inactive. S* indicates that faint activity was later detected despite a stellar appearance. 

^ Note that 8P/Tuttle is not a JFC, but a NIC with a short period; it has an inclination of ~ 55° and Tj = 1.6. 

e 2006 BZ8 was observed to search for activity in this newly discovered object, which has a comet-like (Tj = —1) orbit. None was found. 
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be assumed to have constant intrinsic brightness. Calibra- 
tion of these differential magnitudes onto a st a ndard scale 
was achieved by taking observations o f lLandoltl l|l992h stan- 
dard star fields, and solving for zero points, extinctions and 
colour terms in the usual way. These were used to calculate 
the magnitudes of the field stars in each frame; taking the 
mean of these values gave a very accurate measurement of 
the brightness of the comparison stars, mg(s) ± a R (s). The 
full calibrated comet magnitude is then given by 

m R {c) = [m' R (c) - m' R (s)] + ^rJs) + k' R 8(V - R) (1) 

where the term in brackets is the differential instrumental 
magnitude (comet - star), and 5(V — R) is the difference in 
(V — R) colour between the star and the comet. When the 
colour term k' R is low (as with the NTT observations with 
\k' R \ < 0.06 during both runs) and the comparison stars 
and the comet have similar colours, the final term is negligi- 
ble. Care was taken to select comparison stars with colours 
in the range expected for cometary nuclei to minimise any 
uncertainties in this colour term, although experimentation 
found this method of obtaining calibrated photometry to be 
robust and independent of which comparison stars were cho- 
sen. In the INT data, with Sloan r' and i! filters, the colour 
terms are significant for even slight differences between star 
and comet colours. As can be seen, the dependance on the 
comet colour in the calibrated magnitude causes the deter- 
mination of both the magnitude and the colours to be an 
iterative process. A direct calibration of the comet photom- 
etry was used to provide a first estimate of its colour (with 
large error bars), and therefore calibrated magnitudes and 
from those a revised (V — R) colour. This method typically 
found a stable solution within three or four iterations. 

It is important to note that while the calibration onto 
standard magnitude scales is necessary for comets in dif- 
ferent star fields, and to allow calculation of absolute mag- 
nitudes and colours for the objects, searching for periodic- 
ity can be done just as well using differential light-curves. 
Where possible the differential light-curve is preferred as it 
contains no sources of error beyond the photometric mea- 
surement uncertainties, allowing optimum measurement of 
the relative variation of brightness with time. 



3 RESULTS FROM SNAP-SHOT IMAGING 

'Snap-shot' images of all potential target comets were taken. 
There are three possible outcomes from imaging of comets; 
the comet is either undetected, unresolved or active: Here 
we present the snap-shot data divided into three sections ac- 
cording to these appearance categories. Undetected comets 
are either simply too faint to see, in which case upper limits 
on the size are found, or confused with background sources, 
in which case very little information can be obtained. Un- 
resolved comets are detected and appear as point sources, 
and may be interpreted as being bare nuclei, depending on 
evaluation of the contamination by dust comae. Those which 
appeared to be unresolved, bright and well separated from 
background stars were followed as time-series (see section[4]); 
results on the remainder are presented here. Active comets 
have resolvable coma in individual images; limits are placed 
on the level of activity and the size of the nucleus. All these 
results are summarised in table [2] 



For those comets which were detected, a surface bright- 
ness profile is measured to search for faint coma or quan- 
tify the level of activity. Images and profiles are shown in 
figures [1] and [2] Any coma present is assumed to be in a 
steady state, where the surface brightness is inversely pro- 
portional to the distance from the image centre p (arcsec). 
Fo r apparently inactiv e come ts, we use the equation given 
bv lJewitt fc Danielsor] (|l984 ) 

m c (p) = E c (p) - 2.5 log(27rp 2 ). (2) 

to give the flux from the coma m c within p based on the 
surface brightness E c at this distance. The fraction of the 
flux from any faint coma can then be found by comparing 
m c with the total flux within p. When the fraction is low, 
it can be assumed that the comet is effectively inactive, and 
that the nucleus has been observe d. For active comets we 
make the usual Afp measurement l|A'Hearn et al.ll"l984l ; see 
also paper II), and place limits on the nucleus. 

We obtain the radius of the nucleus tn (or upp er limits 
on it ) based on the observed magnitude m R using l|Russelll 

MB) 

p R r% = 2.238 x 10 22 ^A 2 # 4(me " mRWa) (3) 

or, in terms of the absolute magnitude m R (l, 1, 0), 

p R rl = 2.238 x lo^io - 4 ^©-" 1 ^ 1 ' 1 ' " (4) 

where p R is the geometric albedo, R^ and A the heliocentric 
and geocentric distances in AU, m© = —27.09 is the appar- 
ent _R-band magnitude of the Sun, a is the phase angle in 
degrees and (3 is a linear phase co-efficient in magnitudes 
per degree. Following our previous papers, we assume p R — 
4% and /3=0.035 for these unknown parameters. 

3.1 Undetected comets 

3.1.1 72P/Denning-Fujikawa 

The radius of comet 72P/Denning- Fujikawa was unknown 
before the observations were obtained therefore we could 
not reliably predict its brightness. Five r'-band exposures 
were taken, giving a total exposure time of 250 seconds, 
but the comet was not detected, even in a combined image. 
Historically, this comet has proved elusive: It was first dis- 
covered in October 1881, but lost until a chance recovery 
during a close pass to the Earth in 1978, which was linked 
to the earlier comet despite it being unobserved during any 
of the 10 intervening apparitionfl Subsequent attempts to 
recover the comet during its 1987, 1996 and 2005 appari- 
tions have failed; the n o n-det ection reported here continues 
this trend. lLamv et al. (2004) quote an unpublished radius 
estimate from the updated catalogue of Tancredi et al. of 
0.8 km, but it is far from clear what data this estimate is 
based on. If it were true, then this radius implies a predicted 
m R — 22.4 at the time of the INT run, which should have 
been just detectable. 

2 The history of previous observations of comets are taken from 
Gary Kronk's Cometography (http://cometography.com/). 
Dates of last observation and orbital parameters 
are also sourced from the Minor Planets Center 
(http://cfa-www.harvard.edu/iau/mpc.html) and the JPL 
small bodies database (http://ssd.jpl.nasa.gov/sbdb.cgi ) 
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There are three likely interpretations of this non- 
detection; firstly, the orbit could be incorrect, and in fact 
Denning in 1881 and Fujikawa in 1978 discovered two dif- 
ferent comets. Assuming that this is not the case, the alter- 
native explanations are that the comet has ceased to exist 
(i. e. broken up) since 1978, or that it is a very small nu- 
cleus that exhibits little out-gassing near perihelion, and has 
consequently escaped detection despite numerous attempts. 
If the comet exists, and was in the predicted position, then 
the fact that no object was observed to jur > 23.0 implies 
an absolute magnitude for this comet of mn(l, 1,0) > 18.3 
and therefore tn < 0.6 km. 



3.1.2 75P/Kohoutek and P/2001 X2 (Scotti) 

Both 75P/Kohoutek and P/2001 X2 (Scotti) were observed 
during the 2005 INT run at right ascension 18 h - 19 h 
and declination ~ -20°, which unfortunately placed them 
in the galactic plane [galactic co-ordinates (l,b) ~ (10,0) 
- (20,-10)]. This meant that both were observed against a 
dense star field, making identification of the comets diffi- 
cult. For P/2001 X2 this was exacerbated by the loss of one 
of the frames taken due to a malfunction in the CCD cam- 
era, leaving only two r'-band frames. It proved impossible 
to identify either comet, even with the aid of astrometric 
predictions of their precise positions within these fields. As 
there is a high likelihood of the comet being blended with 
one of the stars within the predicted area it is also impos- 
sible to constrain the brightness of the comets other than 
to say that they must be considerably fainter than the field 
stars, but this does not put meaningful limits on the radius 
as for any sensible size the comet would fit this constraint: 
for 75P the stars in the search area have magnitudes in the 
range van « 14.5 — 20.5, while the comet is exp ected to have 
tur > 22.8 based on a limit of tn < 1.6 km l|Lowrv et al.l 
2003 A from a non-detection in 1999. 



3.1.3 135F '/Shoemaker-Levy 8 

135P/Shoemaker-Levy 8 was inbound at _Rh = 3.6 AU at 
the time of the 2006 INT run, on its way to perihelion at 
2.7 AU in May 2007, although it has yet to be recovered (as 
of 11 July 2007). The comet was not detected in any of the 
five 80 second r'-band exposures or in a combined frame. 
There is no object at the position given by an astrometric 
fit to m_R > 23.5, giving a limit to the absolute magnitude 
of m_fj(l, 1, 0) > 18.4 and therefore a radius of tn < 0.6 km. 



3 Note that the R-band mag nitude of the Sun u se d by 
Lowrv et all dl999l . |2003|) and lLowrv fc Fitzsimmonsl ll200li , 
20051 ). mm = —27.26, w as bas e d on the often quoted value of 
(V - R)q = 0.52 from lAllenl <1973h . which is for the John- 
so n filter system, not the modern Johnson-Cousins system used 
by lLan dolt (199 2]). In the Johnson-Cous ins system (V — R)q = 
0.35 llHolmbere. Flvnn fc Portinaril2006h . The difference between 
ttiq = —27.26 and the true value of m,Q = —27.09 means that 
while the magnitudes measured in these works are accurate, the 
true radii are io°- 2 (- 27 09 + 27 26 ) = 8% larger. The radius for 
75P, and for all others quoted from this group of papers, have 
been corrected to account for this. 



3.14 P/2004 Tl (LINEAR-NEAT) 

P/2004 Tl (LINEAR-NEAT) was discovered in September 
2004 before passing through perihelion at 1.7 AU in Novem- 
ber 2004. It was last observed outbound in March 2005, and 
was still outbound at 3.8 AU during February/March 2006. 
No further information has been recorded on this comet. 
5x60 second r'-band exposures were taken, but the comet 
was not detected. An astrometric fit and measurement of a 
limiting magnitude of mn > 23.6, imply rriR(l, 1,0) > 18.3 
and rN < 0.6 km. 



3.2 Unresolved comets 

3.2.1 SP/Tuttle 

The results from the snap-shot observations of 8P/Tuttle 
obtained during the July 2005 INT run are presented here, 
although they are not included in the analysis of the ensem- 
ble properties of ecliptic comets. 8P is an unusual comet; by 
the 'classical' definition, it is a JFC, as it has an orbital pe- 
riod less than 20 years (P = 13.5 yr), howeve r on dynamica l 
grounds it is a Nearly Isotropic Comet (see iLeviso 3 ll996h . 
as it has a Tisserand parameter with respect to Jupiter of 
Tj = 1.6, and has a Halley-like high inclination (i = 55°) 
orbit. It therefore is likely to be an Oort cloud comet which 
has been captured into a shorter period orbit by interac- 
tions with the planets. It is also known to have a reasonably 
lar ge (rN = 7.8 km) nu cleus, from snap-shot observations 
bv lLicandro et all {2000), and was therefore predicted to be 
a bright target suitable for time-series observation at — 
7.4 AU with the INT, providing an opportunity to add to 
the very few data on the rotational properties of Oort cloud 
comets. However, upon observation the target field was seen 
to be far too crowded for time-series observation, and at- 
tempts to observe the comet during this run were abandoned 
after only a short series of 13 r'-band frames. It was possi- 
ble to separate the comet from background sources in the 
majority of frames; the average ran = 21.61 ± 0.03 is taken 
over the 9 of these frames over which the magnitude of the 
comet was ~ constant, with a brightening in the remaining 
4 being attributed to proximity to a background star. This 
magnitude gives mn(l, 1, 0) = 12.88, corresponding to a ra- 
di us of 7.58±0.12 km, in excellent agreement with the value 
of iLicandro et al.l (|200Ch . 

3.2.2 UP/Remmuth 2 

A snap-shot of 44P/Reinmuth 2 was published in paper II. 
It was observed again during the 2006 INT run, when a total 
of 5 r'-band images were taken over the course of the two 
nights. The comet was detected on both nights, although was 
very faint in poor seeing on the first. The comet's magnitude 
was constant, within the error bars on each point, over the 
two nights at ran — 22.49 ± 0.06, and apparently inactive. 
The profile is consistent with the background PSF, within 
the error bars on each point, and equation [5] gives m c > 
23.26, or < 50±72% of the flux. The absolute magnitude was 
measured to be m,R(l, 1,0) = 15.95, implying rN = 1.84 ± 
0.05 km, which is consistent with our previous measurement 
of r N = 1.96 ± 0.11 km. Colours of (V - R) = 0.31 ± 0.09 
and (R — I) — 0.42 ± 0.13 were measured. 





log(rhn) (nrcsac) Ingfrho) (arcaec) 



Figure 1. Images and surface brightness profiles for unresolved snap-shot targets. From the top left, these are 8P, 70P, fl4P, 120P, 
131P, 160P, P/1995 Af, P/2004 H3, and the Damocloid 2006 BZ8. Notes: 70P - This is a single 295 s r'-band frame; in subsequent 
frames the comet moved towards the two stars and was not separable from them. The profile shows the effect of the nearby star; no 
sensible constraints on the activity level of this comet can be made. 114P - taken on A/7/05, a sum of 8x 90 s r'-band frames. The 
background stars are trailed in this image as it has been corrected for the comet's motion; with only 8 frames it was not possible to 
construct a good background image to subtract the stars. The profile is inactive in the inner part, although dominated by sky noise and 
the star trails in the outer section. 120P - taken on 2/3/06, which is seen to be very faint, but distinguishable from the trailed stars in 
this combination of 4 r'-band frames. The star trails give two additional peaks in the profile at larger p, although the inner 1" implies 
a bare nucleus. 131P - as a faint but star-like source in a combination of 3 frames taken on 2/3/06. The profile is in agreement with 
the image PSF within the error bars, implying that the comet was inactive. 160P — which was detected at a very faint level, even in 
this combined image of three 140 s r'-band frames taken on 2/3/06. The profile appears stellar in the inner part, but is confused with a 
background galaxy in the outer part. P/1995 Al - which was recovered on 2/7/05 returning to the inner solar system 10 years after it 
was last observed. The orbit is accurate, as the predicted motions allowed the 8 frames to be combined to give a star-like comet image 
(with trailed background stars). The profile matches the PSF, even in the outer noisy part, except at p ss 5" where flux from the nearby 
star trail raises the measured surface brightness. P/2004 H3 - taken on the 5th of July 2005, in which the comet is just detectable as 
a very faint point source among the trailed stars. The profile is noisy due to the extreme faintness of the comet, but is consistent with 
an inactive nucleus. 2006 BZ8 - This Damocloid asteroid was observed on 2/3/06 and was bright in individual 10 s exposures, and is 
clearly stellar in this combination of all 15 r'-band exposures. The profile confirms this. 



3.2.3 70P/Kojima 

lLamv et alj (|2004h quote measurements for 70P/Kojima 
from their own unpublished work, from a 'partial rotational 
light-curve', although no indication is given of the quantity 
of data involved. Accordingly, their results (r^ = 1.86 km, 
a/b > 1.1 and P rot > 22 hours) are treated with caution, 
although the radius measurement can be taken at face value 
if regarded as a snap-shot. Identification of the comet at the 
telescope was hampered by the faintness of the target and 
its proximity to a bright field star; only in the first frame 



was it sufficiently separated from the star to allow photom- 
etry. Even in this first frame a small amount of flux from 
the star may have been present within the aperture centred 
on the nucleus; the following radius should be considered 
to be an upper limit. The measured comet magnitude is 
m_R = 22.53±0.15, which corresponds to mn(l, 1, 0) = 15.59 
and rN < 2.18 ± 0.15 km. Under the assumption that this 
is unaffected b y star light, and taken with the radius mea- 
surement from lLamv et all , this im plies a/b > 1 .4, which is 
typical for JFCs. This is larger than lLamv et alj 's value, but 



Observations of 23 distant JFCs 7 



even if they did observe the full range in magnitudes in their 
partial light-curve then a change of viewing angle could in- 
crease the range observed at a different epoch: light-curves 
only give a minimum a/b when the orientation of the rota- 
tion axis is unknown. Based on both data sets, an average 
tn = 2.0 ± 0.2 seems a reasonable value to assume for the 
size of 70P. 



3.2.4 n 4 p / 'Wiseman-Skiff 

114P/Wiseman-Skiff was observed with the INT in July 
2005 when inbound at 3.75 AU, and was detected but at 
too faint a level to produce a reliable light-curve. The com- 
bined image of all 8 frames and corresponding profile show 
that the comet appeared inactive. The profile was only mea- 
surable to p — 4", but gives m c — 24.45 and therefore a 
limit on the coma contribution of < 25 ± 41%, with the 
large error bars due to the noisy profile in the outer part. 
The measured m R = 22.95 ± 0.11 gives m R (l, 1, 0) = 17.34 
and tn = 0.97 ± 0.05 km. This agrees with the as sessment 
that t he nucleus is a relatively small one from lLamv et al.l 
(2004), who again quote their own unpublished result of 
tn = 0.78±0.04, based on an HST snap-shot when the comet 
was at 7?h =1.6 AU and presumably active. Assuming that 
these two snap-shots were taken at opposing light-curve ex- 
trema, then the minimum required axial-ratio is a/b > 1.5, 
and the average effective radius is tn = 0.89 km. 



3.2.5 120P /Mueller 1 

Although not detectable in individual images, 
120P/Mueller 1 was found in a co-added image of 4x 
115 s r'-band frames shifted to account for its predicted 
motion. The surface brightness profile for this comet 
matches the stellar PSF within the inner arcsecond, al- 
though beyond this there are two further peaks due to the 
nearby trails of background stars, making the formal limit of 
m c > 23.1 meaningless. The comet had tur — 23.52 ± 0.13 
in the combined image. This gives m,R,(l, 1,0) = 17.84 and 
r\- = 0.77 ± 0.05 km. Assuming that this measurement is 
of the radius of the bare nucle us, it can be compar ed with 
the previous measurement of lLowrv et al.l (| 19991 ). They 
measured tn = 1.7 ± 0.2 with a snap-shot of a star-like 
nucleus at = 3.1 AU, which at approximately twice the 
radius measured here is a large discrepancy; such disparate 
measurements would require a/b = 4.7 to be explained 
purely by the shape of the nucleus, which is rather high. It is 
p ossible that th ere was either a non-zero coma contribution 
in lLowrv et alj 's -Rh = 3.1 AU measurement which was not 
present at 3.9 AU, or that the comet has a steeper tha n 
normal phase function / opposition surge l|Lowrv et alj 's 
observations were taken at a — 2.9°). The unknown colour 
of the nucleus could also contribute, but only at a level 
similar to the uncertainty on tur: If 120P has a very re d 
nucleus the difference in colour term between lLowrv et alj 's 
i?-band and our r'-band measurements could cause the tur 
presented here to be up to a tenth of a magnitude too faint. 
It is likely that some combination of these factors, combined 
with a more modest elongation of the nucleus, is required 
to explain the differences between these measurements. 



3.2.6 13 IP/ Mueller 2 

131P/Mueller 2 was outbound at 3.5 AU when observed with 
the INT in 2006. Three r'-band frames were taken, and the 
comet was found and observed to have a faint, star-like ap- 
pearance. The large error bars on the surface brightness at 
large p make the coma limit of m c > 23.1 (< 68 ± 94% of 
the total flux) almost meaningless, but the profile appears 
star-like and the comet can be assumed to have been inac- 
tive. The average magnitude was tur = 22.77 ± 0.09, corre- 
sponding to an absolute magnitude of m.R(l, 1, 0) = 17.59, 
and therefore a radius of rN = 0.87 ± 0.04 km. In addition, 
three U-band frames were taken, allowing measurement of 
the colour of the nucleus. A value of (V — R) = 0.45 ± 0.12 
was obtained, which falls in the centre of the observed JFC 
distribution (see section [6}. 



3.2.7 160P /LINEAR 

Observations of 160P/LINEAR from the INT during 2006 
detected the comet when outbound at ~ 4 AU. The comet 
appeared as a very faint but apparently inactive object in 
3x 140 s r'-band exposures. The presence of a nearby galaxy 
prevents measurement of limits on any coma, with m c (5") = 
22.6. The inner part of the profile is a good match for the 
background PSF; measuring the surface brightness at 1.3" 
gives m c > 28.0, < 1.9±5.8% of the total flux. The measured 
value of m,R — 23.69 ± 0.18 implies mfl(l, 1, 0) = 17.58 and 
r N = 0.87 ± 0.07 km. 



3.2.8 P/1995 Al (Jedicke) 

An attempt to detect the return of P/1995 Al (Jedicke) 
was made in July 2005 with the INT, although the comet 
was still at Rh = 5.5 AU. A faint object was discovered 
near the nominal position with an average rriR = 22.63 ± 
0.07. Measured in a combination of all 8 r'-band frames, 
the comet's surface brightness profile appeared stellar in the 
inner region, although was dominated by sky noise beyond 
2", which is taken to show that it was inactive, as would be 
expected with the comet still beyond Jupiter's orbit. This 
implies a radius for the nucleus of rN = 2.48 ± 0.08, and an 
absolute magnitude of mfl(l, 1,0) = 15.31. 



3.2.9 P/2004 H3 (Larsen) 

P/2004 H3 (Larsen) was only at its second opposition when 
observed with the INT in July 2005, having been discovered 
the previous April. It was at R^ — 3.7 AU, outbound, and 
was only just recovered in a frame combining all six 90 s r'- 
band exposures. The measured magnitude, in the combined 
frame, was rriR — 24.30 ± 0.21, implying that the nucleus is 
small, with m_fj(l, 1, 0) = 18.67 and a corresponding radius 
of rN = 0.53±0.05 km. A profile is just measurable, but only 
within 3", and is very noisy. Formally, the coma magnitude 
from equation [2] is m c > 24.6, implying a flux contribution 
of < 80%, although the error on this is also ~ 80%, making 
it of limited value. The profile matches the stellar PSF at 
all radii, within the error bars. 
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Figure 2. Images and surface brightness profiles for 56P, 78P and 
103P. The comets are visibly active. The solid diagonal lines in 
the top right corner of each profile show gradients of -1 and -1.5, 
from two theoretical models of steady state coma; the comets all 
have profiles with gradients ~ —1.5. 



3.2.10 2006 BZ8 

Observations of 2006 BZ8 were carried out to search for ac- 
tivity around this newly discovered object. It was regarded 
as being almost certainly a comet, with Tj — —1.0. No 
activity was observed though, and the object was bright 
(m,R = 16. 858± 0.003) and stellar. A surface brightness pro- 
file entirely matched the combined frame PSF, and equation 
[2] gave a limit on any coma at m c > 21.67, or < 1.2 ± 1.2% 
of the flux. 2006 BZ8 is officially a CentauiQ as it has an or- 
bital semi-major axis of a = 9.7 AU, however it seems likely 
that it is from the Oort cloud due to its Tisserand parameter 
and high inclination (i — 165°), an d is better cl assified as 
a member of the Damocloid family (|jewittl [20051 ). The ab- 
solute magnitude was measured to be mn(l, 1,0) = 14.14, 
corresponding to tn = 4.22 ± 0.02 km, and its colours were 
found to be (V - R) = 0.62±0.01 and (i?-/) = 0.36±0.01. 
These results are not included in the discussion on the en- 
semble properties of JFCs in section [6] 



3.3 Active comets 

3.3.1 56P/ Slaughter-Burnham 

56P/Slaughter-Burnham was discovered in 1 958, and has 
previ ously been found to have t-n = 1.56 km (M eech et al.l 
|2004| ). This led to a predicted nuclear magnitude of tur = 
22.1 at the time of the 2006 INT observations, and the comet 



4 Objects with 5.5 < a < 30.1 AU are designated Centaurs. 



was well placed for time-series observations at the begin- 
ning of the night. The comet was active at Rh — 3.8 AU 
(outbound). The mean magnitude was txir = 20.73 ± 0.04, 
implying m_fj(l, 1,0) > 14.90 and tn < 3.00 km, i n broa d 
agreement with the measurement by I Meech et al.l l|2004h . 
The colours of the comet were measured to be (V — R) = 
0.51 ± 0.05 and (R - I) = 0.35 ± 0.08, which fall in the 
observed range for nuclei, although it should be noted that 
these colours represent the colours of the coma, as the comet 
was active. These colours are similar to ot her measure- 
ments of th e broad-band colours of dust comae |Lowrv et alj 
1 19991 . 120031 ), and also appear to be in agreement with the 
trend of cometary dust to have more neutral (less red) 
colours at longer wavelengths, which is interpreted as show- 
ing that the dominant grains are larger than 1 /im in size 
l|jewitt fc M eech 1983). It cannot be assumed that the coma 
has the same colour as the surface of the nucleus for highly 
active comets; even if the dust component is of similar com- 
position, there is likely to be gas present also, especially in 
the U-band. The flux within an aperture of p — 5"(= 11400 
km at the comet) gives Afp = 22.3 ± 0.3 cm. 

3.3.2 78P/Gehrels 2 

Three r'-band frames were taken of 78P/Gehrels 2 during 
the 2006 INT run, which showed the comet to be clearly 
active when outbound at Rh — 3.8 AU. lLowrv fc Weissmanl 
(2003) measured a radius of = 1.54±0.12 km from a snap- 
shot taken when the comet was at Rh — 5.5 AU; the upper 
limit obtained here of < 4.21 km, from ran = 19.42±0.01 
and mn(l, 1, 0) > 14.16, is consistent with this. The activity 
level, measured within p — 5"(= 10400 km), was found to 
be Afp= 18.2 ±0.2 cm. 

3.3.3 103P/ Hartley 2 

103P/Hartley 2 was observed to be highly active during the 
2005 NTT run (see paper II). It was observed again dur- 
ing the 2006 INT run, despite a strong suspicion that it 
may still be active at Rh = 5.0 AU, as at the time this 
comet was a potential target for the Deep Impact extended 
mission, The comet was detected in two r'-band 

frames taken on the first night, and was clearly active. The 
measured tur — 21.00 ± 0.05 implies a radius upper limit 
of tn < 4.45 km (it is thought that the nucleus is actually 
much smaller; see paper II), while the flux within p — 15600 
km gives Afp — 47.2 ± 0.7 cm. The coma had colours of 
(V-R) = 0.16 ±0.09 and (#- I) = 0.35 ±0.08. The rather 
blue (V — R) measurement could reflect the presence of sig- 
nificant flux from gas emissions within the U-band, such as 
the C2 bands at 5165 and 5635 A, while the (R — I) colour 
is similar to that found for nuclei and dust comae. 



4 TIME-SERIES PHOTOMETRY 

Those comets which appeared (upon visual inspection at the 
telescope) to be inactive and reasonably bright were followed 

http://epoxi.umd.edu/ - 103P has now been confirmed as the 
target of this mission, with a fly-by date in 2010. 
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Table 2. Derived values and limits on radii and activity from snap-shot photometry. 



Comet van mfl(l,l,0) fN [km] Afp [cm£ 
UNDETECTED 

72P > 23.0 - > 18.3 < 0.6 < 0.45 
75P _____ 

135P > 23.5 - > 18.4 < 0.6 < 0.88 

P/2001 X2 - - - 

P/2004 Tl > 23.6 - > 18.3 < 0.6 < 0.57 
UNRESOLVED 



8P 


21 


.61±0. 


03 


>25.1 (4.3") 


12 


.88±0. 


03 


7. 


58±0. 


12 




44P 


22 


.49±0. 


06 


>23.3 


15 


.95±0. 


06 


1. 


84±0. 


05 


<1.7 


70P 


22 


.53±0. 


15 


>23.9 


15 


.59±0. 


16 


2. 


18±0. 


15 


<23.7 


114P 


22 


.95±0. 


11 


>24.5 (4.0") 


17 


.34±0. 


11 


0. 


97±0. 


05 


<0.24 


120P 


23 


.52±0. 


13 


>23.1 


17 


.84±0. 


13 


0. 


77±0. 


05 


<2.9 


131P 


22 


,77±0. 


.10 


>23.1 


17 


.49±0. 


10 


0. 


87±0. 


.04 


<0.50 


160P 


23 


.69±0. 


18 


>28.0 (1.3") 


17 


.58±0. 


18 


0. 


87±0. 


07 


<4.7 


P/1995 Al 


22 


.63±0. 


07 


>22.1 


15 


.31±0. 


07 


2. 


48±0. 


08 


<5.7 


P/2004 H3 


24 


19±0. 


21 


>24.6 (3.0") 


18 


.55±0. 


21 


0. 


56±0. 


06 


<0.69 


2006 BZ8 


16 


.88±0. 


01 


>21.7 


11 


.16±0. 


01 


4. 


22±0. 


.02 


<41.4 



ACTIVE 

56P 20.73±0.04 - >14.90 <3.00 22.2±0.3 

78P 19.42±0.01 - >14.16 <4.21 18.2±0.2 

103P 21.00±0.05 - >14.04 <4.45 47.2±0.7 



a Coma magnitude within p = 5" (unless stated otherwise) measured using equation [2] 
^ Afp measured through an aperture of radius 5". 



over an extended period in order to obtain a light-curve. We 
assume that the variation in brightness is due to the nucleus 
being a rotating non-spherical body, and that the variation 
in reflected flux is due to the varying projected surface area. 
Here we follow the common practice of describing the nu- 
cleus as the simplest non-spherical shape: a tri-axial ellip- 
soid with semi-axes a > b = c. We therefore expect to see a 
double peaked light-curve with a rotation period P ro t twice 
the fitted Fourier period Patted- The Fourier method used 
to search for periodicities is described in paper I; briefly, 
we search for minima in the range x 2 l v = 1 ± v/2/i/ in a 
reduced \ 2 periodogram. v is the number of degrees of free- 
dom of the model, given by v — (N — 3) where N is the 
number of data points, for a first order Fourier model. 

The periodicity search allows us to constrain the rota- 
tion period, while the peak-to-trough amplitude of the light- 
curve Am gives the minimum elongation of the nucleus via 
a/b > 10 0,4Am . We can only measure a lower limit on a/b, 
as we do not know the orientat ion of the rotation axis . We 
use the approximation given bv lPravec fc Harris! (|2000T ): 



to calculate the minimum bulk density of each nucleus us- 
ing the light-curve derived parameters a/b and P ro t (hours). 
Note that this density value is a minimum for two reasons: it 
relies on a/b, for which we have only a lower limit, and also 
that the nucleus need not be spinning at its break up rate, 
but clearly cannot be spinning faster than it. The results 
from the time-series observations are summarised in table [3] 

4.1 40P/Vaisala 1 

40P/Vaisa,la 1 was observed over 7 nights during the 2005 
INT run. It had a mean m R = 22.08 ± 0.02. The sur- 
face brightness profile is inconclusive, but we suspect that 
the comet was weakly active at = 4.6 AU. The coma 
contribution, assuming steady state, was measured to be 
m c > 22.22 within 5", and therefore up to ~ 95% of the 
flux, but the irregular shape of the profile (fig. [3} implies 
that the steady state assumption may not be valid in this 
case. 

With observations over 7 nights, changes in the relative 
positions of the Earth and the comet meant that there was 
a non-negligible change in observing geometry. The change 
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Figure 3. Image showing 40P, made up of 20x85 second expo- 
sures taken on the 6th night of the 2005 INT run. The profile is 
inconclusive, but we suspect that the comet was weakly active. 
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Figure 4. Pcriodogram for 40P. No convincing period was found. 

in apparent magnitude due to variations in th, A and a be- 
tween the first and last frames is Smn ~ 0.04 mag., which 
although small compared with the uncertainty on individ- 
ual points (typically ~ 0.2 mag.) is considerable when com- 
pared to the error on the average magnitude. The photom- 
etry was therefore reduced to the appropriate ma (1,1,0) 
using the precise position of the comet at the time of each 
observation, which was generated using HORIZON^. The 
phase function is assumed to be linear with j3 = 0.035 
mag. deg -1 , as the range in phase angle over 7 nights is 
not large enough to independently measure (3 for 40P. The 
mean mn(l, 1, 0) = 15.72 gives an upper limit to the ra- 
dius of tn < 2.05 ± 0.02 km, assuming a 4% albedo. This 
is a stronger cons traint than the pre vious upper limit of 
tn < 3.6 km from lLowrv et all |l999), who did not detect 
the comet when it was at Rh = 6.01 AU, outbound, in 1995. 
If the coma contamination is large, as suggested above, then 
the true radius of the nucleus could be considerably smaller 
than this. 

The absolute magnitude variations were searched for 
6 http : //ssd. j pi .nasa.gov/horizons . cgi 



47P/AshbmQk-Jacksan [R] 




Figure 5. Images of 47P taken during the 2006 INT run. There 
is a faint tail visible to the West (left in these images) of the 
comet, which is clearer in the night 3 data (top). The profile shows 
activity. In the night 4 data (bottom), the tail is less obvious, as 
the total exposure time is less. 

periodicities, but none were found which were convincing. 
The strongest minimum in the periodogram (fig. |4} corre- 
sponds to Petted = 15.3 hours, and therefore a long rotation 
period of 30.6 hours, but does not give a visually acceptable 
light-curve at either of these periods. Periodogram searches 
were also carried out on short subsets of the data, using 
fully differential light-curves for each pair of nights, during 
which the comet's motion was small enough to give common 
comparison stars. Shorter periods in these data are neither 
particularly convincing in the subsets, due to sparse data, 
nor good fits to the full data set. Although there is a large 
range in observed magnitudes (Am ~ 1.4 ± 0.3 mag), this 
appears to be due to variations in the dust coma, and is 
not due to the nucleus, again implying that in this case the 
coma completely dominates. 

The colour measurements have large uncertainties due 
to the comet being very faint, but are consistent with means 
of (V — R) = 0.37 ± 0.10 and (R - I) = 0.65 ± 0.11. At 
large distance from the Sun the coma should be dominated 
by dust, with very little gas present; the colours measured 
for the comet could therefore be similar to the colours of 
the nucleus, although this is actually a measurement of the 
colour of the dust coma. 

4.2 47P/Ashbrook- Jackson 

A time series on 47P/ Ashbrook- Jackson was published in 
paper II. The comet was observed again during the 2006 INT 
run; a total of 5 r'-band exposures totalling 18.5 minutes 
were taken over the two nights. Upon processing these data, 
47P was found to be faintly active, despite still being beyond 
5 AU, with a faint tail visible to the West of the comet (fig. 
[5]). The activity is clearly weak: The profile measured in 
good (~ 1") seeing on the last night matches the stellar 
PSF within the inner part, where the nucleus appears to 
dominate the flux. The activity level can be quantified as 
Afp = 3 — 11 cm, with the lower number being measured 
on the latter night. 

The average magnitude over the two nights is rriR — 
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Figure 6. 47P mean magnitudes from paper II, the INT data 
and from the literature reduced to rjj = A = 1 AU and plotted 
as a function of phase angle a. There is a good fit through the 
time-series data giving a linear phase law with (3 = 0.083. Labels 
give the sources o f each point; LiOO: Licandro ct al. 2000); La04: 
lLamv et a"fll2004l ; NTT: paper II. The faint activity seen in the 
INT data means that it provides a lower limit on m R (l,l,a), 
although the coma clearly contributes little flux as making this 
point any fainter would mean an even steeper phase function and 
a poor fit. 




Figure 7. Co-added image of 24 frames taken of 94P on July 
6th, 2005. Each frame had an exposure time of 75s, giving this 
combined frame an equivalent exposure time of ~30 minutes. The 
profile implies that the comet could have been weakly active at 
the time of observation, although the rise at the end is more likely 
due to residual star light. 



data for 47P (including that of iLicandro et alb gives (V 
R) = 0.42 ± 0.02, (R — I) = 0.44 ± 0.03. 



22.01 ± 0.04, with a total variation similar to that observed 
in the light-curves above, of Am = 0.4 mag. These imply 
m R (l, 1, 0) > 14.93 and r N < 2.96±0.05 km, with a/6 > 1.4, 
slightly fainter than our earlier results despite the weak 
coma, possibly due to a steeper phase function than /3 — 
0.035 mag. d eg" 1 . Taking bo t h the NTT and INT resu lts, 
and those of ILicandro et alj (l2000l ) and lLamv et all 7 ! , al- 
lowed calculation of a rough phase function (fig. [6}. A best fit 
straight line gives a steep phase function of /3 = 0.083±0.006 
mag. deg -1 , implying a true m R {l, 1,0) = 14.477 ± 0.007 
and consequently rN = 3.63 ± 0.01 km. The fit to the three 
time-series measurements (taking the 5 points over 2 nights 
from the INT data as a 'time-s eries') is excellent, while the 
snap-shot from ILicandro et al.l also agrees with this fit if it 
is assumed that it was taken near a light-curve minimum. 
Therefore we regard this as a reasonable estimate of the 
phase function of 47P, despite the fact that it is unusually 
steep. 

Colours were measured of (V — R) = 0.29 ± 0.06 and 
(R — I) = 0.79 ± 0.07. These are unusual when compared to 
typical nuclei colours, an d show a large change from those 
presented in paper II and ILicandro et all s (V — R) = 0.4 ± 
0.3. 47P was weakly active at the time of these observations, 
so we favour our earlier value. A weighted mean of all colour 



In this case taking the rN=2.8 km quoted bv lLamv etafl (2004) 
and calculating the nuclear magnitude at the time of their HST 
observations. The uncertainty on this radius is not reported: The 
error bar on this point is set equal to that on INT data (0.04 
mag.) for the purposes of weighting the best fit. This does lead 
to some uncertainty on the phase function; if the error bar is 
instead set equal to that from our extended NTT observations 
(0.007 mag.) then the corresponding increase in weight on the 
point gives /3 = 0.076 mag. deg -1 . 



4.3 94P/Russell 4 

94P/Russell 4 was observed over 4 nights in July 2005 using 
the INT. A total of 83 r'-band frames were taken, and in 
addition at least one colour block on each night. The com- 
bined image and corresponding profile are inconclusive (fig. 
O; the profile shows some signs of faint activity beyond 2" 
from the nucleus, although this may be due to residual flux 
in the combined image from particularly bright nearby stars. 
The calculated coma contribution within 5" is m c > 21.25, 
or < 77±30%; measured before the peak from possible stel- 
lar contamination the contribution is m c (3") > 22.06, or 
< 37 ± 15%. 

The mean magnitude of 94P was measured to be m R — 
20.974 ± 0.015; as with other comets observed during this 
run the individual magnitudes were converted to absolute 
magnitudes due to the 5m R ~ 0.03 mag. difference over 
the four nights, primarily due to the changing phase an- 
gle. The absolute magnitude and radius were measured to 
be m R (l, 1,0) = 15.187 and r N = 2.62 ± 0.02 km. These 
are technically upper limits if there was any weak activity, 
however we believe that the comet was effectively inactive 
as these are in excellent agreement with preliminary results 
from a light-curve taken by the authors for a different pro- 
gram when the comet was clearly inactive at 4.7 AU in July 
2007. 

The motion of this comet was such that it was not pos- 
sible to use common stars to produce multi-night differential 
light-curves, but since it was reasonably bright and exhibited 
a large amplitude variation the increase in size of the error 
bars due to calibration onto the Landolt scale was not signif- 
icant. The null hypothesis of constant brightness is rejected 
at a 50a level. The periodogram for 94P is shown in fig. [H] 
the strongest minimum is at Pfiued = 16.7±0.4, correspond- 
ing to a rotation period of P Iot = 33.4T0.8 hours. Although 
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Figure 8. Pcriodogram for 94P. 
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Figure 9. Folded light-curve for 94P, period = 33.4 hours. Data 
from separate nights are plotted with different symbols (circle 
- night 1; cross - 2; square - 3; triangle - 4): the variation in 
brightness is considerable even during single nights. 

this minimum only has x 2 jv = 193/79 = 2.45 = 9<r, and 
the data folded onto this period displays considerable scat- 
ter (fig. O, it is clear that there is substantial variation in 
the brightness of the comet, again implying that the nucleus 
dominates the flux. 

The full range observed in the data is Am = 1.2 ± 0.2 
mag., implying a/b > 3. This is quite an extreme elonga- 
tion, in excess of any other measured nucleus, and exceeded 
by few asteroids. The fact that there are large ranges in ob- 
served brightness on individual nights and that this comet 
was bright show that this is not due to any calibration issue. 

An alternative explanation is that the range in ran is 
due to variations in albedo across the nucleus surface instead 
of changing cross-sectional area, which would be revealed by 
periodicities in the colour indices. There are some variations 
in the measured colours of the comet, although the individ- 
ual measurements are normally distributed around the mean 
values of (V - R) = 0.62 ± 0.05 and [R-I) = 0.44 ± 0.06. 
The variations do not appear to be correlated with either of 
the above periods, nor can the r'-band data be fit to peri- 
odicities found in the (R — I) sequence. Changes in surface 
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Figure 10. Co-added image of 121P, from the INT frames taken 
on 2nd March 2006. The nucleus appears star-like, and the profile 
implies that the comet was inactive at the time of observation, 
yet there is a tail visible ~ West [upwards in this image] of the 
nucleus (highlighted in fig, lilt . 

colour do not seem to be the cause of the large amplitude 
variations seen in the light-curve. 

4.4 121P/Shoemaker-Holt 2 

121P/Shoemaker-Holt 2 was predicted to have m_R = 22.18 
during the 2006 INT run, based on a radius of tn = 
1.75 ±0.63 km from a snap-shot observation in June 1999 by 
lLowrv et al] (|2003l ). The comet was brighter than expected, 
but a combined image and star-like surface brightness profile 
showed the comet to be inactive (fig. 1 10 p , with a calculated 
coma contribution within 5" of m c > 24.95, or < 2 ± 4% 
of the total flux. However, inspection of the co-added frame 
revealed a dust feature in the plane of the comet's orbit. 
Originally suspected to be a dust trail of remnant dust par- 
ticles along the orbit due to the apparent inactivity of the 
comet, it is clear from fig. [11] (in which colours are used to 
highlight faint features, and the orbit of the comet projected 
onto the sky is over-plotted) that the feature is a tail, as it is 
not entirely in the orbital plane and spreads around the anti- 
solar direction. This means that there are two contradictory 
measurements; the stellar nucleus profile implies that there 
is no coma, and there is clearly none resolvable, yet there is 
a tail, which implies that the comet is quite highly active. A 
tempting interpretation is that the observations happened 
to catch the comet just after it 'turned off', or shortly after 
an outburst from an otherwise inactive comet, meaning that 
there is no near nucleus coma, but the material in the tail 
has not yet dispersed. 

Assuming that the nucleus was inactive, the measured 
m R = 20.775 ± 0.006 implies m fl (l,l,0) = 14.660 and 
tn = 3.35 ± 0.01 km. This is co nsiderably larger th an the 
tn = 1.75 ± 0.63 result found bv lLowrv et alj (|2003h . even 
considering the large error bar on their result. This fur- 
ther suggests that there is some activity present in the INT 
frames; the axial-rat io required to explain the difference as 
due to lLowrv et al.l observing an tn = 3.35 km body at a 
light-curve minimum is a/b = 6.3, which is unrealistic. How- 
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(a) (b) (c) 

Figure 11. Wider field co-added images of 121P, using false colour to highlight the faint tail, (a) shows the 1st March 2006 INT data, 
during which the seeing was poor and exposure times were doubled, and there were too few images to generate a background image 
without the comet to subtract. This image is therefore a direct sum of the 4 r'-band images (equivalent exposure time t e = 27 minutes), 
and contains stars, (b) was produced from the 18 r'-band frames (t e = 65 minutes) taken at the INT on the 2nd of March 2006, and has 
had background objects removed, (c) shows the data obtained with the FTN in May 2006, and shows no tail, although it is not such a 
deep image (14 minutes). 



ever, the brightness of the comet does show variations which 
may be periodic and could be interpreted as being due to 
the rotating nucleus: the periodogram and best fit folded 
differential light-curve are shown in fig. [12] and fig. 1131 With 
only 4 frames taken on the first night and 19 on the sec- 
ond, there is ambiguity in the determined period. A number 
of periods give visually acceptable light-curves; Figure 1131 
shows 4 that illustrate the range seen in the periodogram. 
Note that the deepest minima have x 2 / v < 0.68(s -la), 
implying that the error bars on the photometry are over es- 
timated (i. e. the brightness measurements are more accu- 
rate than the conservative estimates on the uncertainties). 
The null hypothesis of constant brightness is rejected, at 
only x 2 /" = 47.4/21 = 2.2 = 3.8a, although this would 
be rejected at a higher confidence level if the error bars 
on individual points were reduced. The total variation is 
Am = 0.15 ±0.03 mag., implying a minimum a/b > 1.1 and 
a fairly spherical shape or near pole-on orientation, although 
if there is near nucleus unresolved coma then the true axial- 
ratio will be larger (see discussion on 36P in section [5}. The 
minimum density implied by this 0.12 ± 0.04 g cm -3 . 

To determine whether or not 121P was still active, and 
in particular to investigate the interpretation that it had 
just ceased outgassing in March, the comet was observed 
again on the 31st of May 2006, using one hour on the FTN. 
Based on the radius measured using the INT, the predicted 
magnitude at this time was ran — 20.69; actually brighter 
than during the INT observations due to the lower phase an- 
gle, and theoretically well within the capabilities of the 2.0m 
FTN. 14 iZ-band frames were taken over the course of the 
hour, along with 3 V-band and suitable standard star obser- 
vations for calibration. As the FTN does not have an auto- 
guider, exposure times were limited to 60 s to minimise trail- 
ing of the images due to the telescope's movement. Combin- 
ing all frames gave a detection of the comet, which was very 
close to the predicted brightness, with ma = 20.63 ± 0.10. 
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Figure 12. Reduced \ 2 periodogram for 121P. 

The profile of this combined image indicated that the comet 
was inactive (m c (3.1") > 24.2, implying a flux contribution 
of < 4 ± 7%), and no tail was detected to a 3a limiting sur- 
face brightness of E.r > 25. Unfortunately this does not rule 
out a tail, as the surface brightness in the centre of the tail in 
the INT data is beyond these detection limits at E_r ~ 25.5. 
The FTN observations were therefore unable to constrain 
the evolution / dissipation of the tail, although the star-like 
profile suggesting inactivity adds support to the idea that 
both sets of observations were of a bare nucleus which had 
recently ceased out-gassing. 

If this is the case then an approximate measurement 
of the phase function can be made. As the measured FTN 
brightness was close to the predicted value, then f3 must 
be close to the canonical 0.035; in fact a fit to these two 
points gives (3 = 0.047 ± 0.020, and m R (l, 1, 0) = 14.494 ± 
0.006, implying tn = 3.61 ± 0.01 km. T h ese re sults are still 
inconsistent with those of lLowrv et al] (|2003t ), and a non- 
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Figure 13. Folded differential light-curves for 121P, periods A : 
8.4, B = 10.1, C = 12.8 and D = 17.4 hours. 



standard phase law does little to improve matters. Including 
their result gives a phase law with /3 = 0.036 ± 0.006, and 
an implied mn{l, 1, 0) = 14.65 and tn = 3.36, but still does 
not fit their magnitude well. The colours measured for 121P 
were (V — R) = 0.53 ± 0.03 and (R - I) = 0.44 ± 0.03, 
from the INT data, and (V - R) = 0.29 ± 0.20 in the FTN 
data, where the extreme faintness in the few short V-band 
exposures gives large error bars. 



4.5 P/2004 H2 (Larsen) 

P/2004 H2 (Larsen) was observed in July 2005 at only its 
second opposition since discovery in April 2004. A total of 
98 r'-band frames were taken; the combined image shown in 
fig. ll4f a1 is made up of the 16 taken on the 4th night. It can 
be seen in this image that the comet was weakly active, with 
a small tail extending to the West of the nucleus. The profile 
also shows the comet to be active, with m c = 23.85 measured 
within p = 5". The comet flux is dominated by the nucleus 
though, as this corresponds to only 13 ± 11% of the total 
flux (Afp = 6.0±0.1 cm, within p = 11, 000 km), under the 
assumption of steady state coma. This assumption is only 
an approximation in this case, as the slope of the profile has 
a gradient slightly steeper than -1.5. 

P/2004 H2 was moving relatively slowly on the sky at 
the time of observation, allowing for both long exposure 
times (210 s) and the use of common stars to generate en- 
tirely differential light-curves over three night time-spans. 
These were searched for periodicities, as was the 7 night cali- 
brated light-curve, adjusted to take into account the chang- 
ing observational geometry. No convincing periodicity was 
identified. 

As the comet was active, we find a lower limit to the 
absolute magnitude and an upper limit to the radius. The 
light-curve does give an accurate measurement of ran = 
21.590 ± 0.007, which reduces to an absolute magnitude of 
?tir(1, 1,0) > 15.875 and therefore a mean effective radius 
of tn < 1.91 ± 0.01 km. The mean measured colours of this 
comet are (V-R) = 0.49±0.02 and (R-I) = 0.28±0.04, al- 
though the individual measurements have considerable scat- 



Figure 14. Image showing P/2004 H2, in which a faint tail can 
be seen to the West (left in this image) of the nucleus. The profile 
also shows that the comet was weakly active, although the flux 
from the nucleus dominates in the photometry. 



ter around these means, possibly due to variations in the 
near nucleus coma. 



5 36P /WHIPPLE 

Time-series observations were taken on 36P/ Whipple on 
three occasions; over 1.5 nights during the March 2005 NTT 
run, over a further two nights a year later with the INT, and 
a final time with the NTT in February 2007. The reason for 
returning to the same comet was that, as with the others ob- 
served during the weather affected NTT run (see paper II), 
its rotation period was not confirmed at the first attempt. 
What made 36P particularly interesting was the fact that 
the rotation period suggested by the first data set was very 
short, at around 3.5 hours. Such a period implied a minimum 
density considerably higher than any other previously ob- 
served nucleus; equation [5] gave a minimum density required 
for gravitational cohesion of 1.2 g cm -3 , greater than that of 
water ice. This pr eliminary result was presented at the 200 5 
ACM conference dSnodgrass. Fitzsimmons fc Lowry||2005al ). 
and was noted bv lToth fc Lissej l|2006). The additional data 
sets do not support this result though, and show 36P's rota- 
tion period to be considerably longer. Here we first describe 
the results from each run before discussing the interpreta- 
tion of the whole data set. 



5.1 NTT data: March 2005 

36P was found to be bright (m_R = 21.37 ± 0.02) and ap- 
peared stellar in each individual frame. A combined image 
(fig. I15p also appears fairly stellar, but there is a slight ex- 
tension towards the East (right in fig. I15f a)) and a sur- 
face brightness profile (fig-HHb)) shows clear signs of coma. 
Equation[5]gives a coma magnitude of m c > 22.4, which cor- 
responds to 40 ± 23% of the total flux being due to steady 
state coma. The flux within an aperture of p = 5" = 12000 
km gives Afp = 20.2 ± 0.3, which is more properly consid- 
ered an upper limit since the slope of the profile is steeper 
than the -1.5 limit for a steady state coma. Such activity at 
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Table 3. Derived physical parameters and colours from time-series photometry on JFCs. 



Comet 


m R 




m«(l,l,0) 


rn 


Prot 






(V-R) 


(R-I) 










[km] 


[hr] 




[g cm" 3 ] 






36P^£ 


21.370±0.016 


>22.4 


15.272±0.016 


2.52±0.02 


3.56±0.02 


1.4±0.1 


1.24±0.07 


0.46±0.05 


0.52±0.04 


36P^ 


21.570±0.008 


>24.5 (3.3") 


15.247±0.008 


2.55±0.01 


~40 


1.9±0.1 


O.OlzbO.Ol 


0.48±0.03 


0.56±0.03 


40F§ 


22.08dz0.02 


>22.2 


> 15.72 


< 2.05 


n/a 


n/a 


n/a 


0.37±0.10 


0.65±0.11 


47F§ 


22.01±0.04 


>22.9 


14.93±0.04 


2.96±0.05 


n/a 


1.4±0.1 


n/a 


0.29±0.06 


0.79±0.07 


94P 


20.974±0.015 


>22.1 (3.0") 


15.19±0.015 


2.62±0.02 


~33 


3.0±0.5 


0.03±0.01 


0.62±0.05 


0.44±0.06 


121F@ 


20.775±0.006 


>25.0 


14.660±0.006 


3.35±0.01 


io±i 


1.15±0.03 


0.12±0.04 


0.53±0.03 


0.44±0.03 


121P£ 


20.63±0.10 


>24.2 (3.1") 


14.60±0.10 


3.44±0.16 


n/a 


n/a 


n/a 


0.29±0.20 


n/a 


P/2004 H^] 


21.590±0.007 


>23.9 


> 15.88 


< 1.91 


n/a 


n/a 


n/a 


0.49±0.03 


0.28±0.04 



a Limiting coma magnitude measured within 5", unless stated. 
^ Lower limits as the orientation of the rotation axis is unknown. 

c Results from 2005 NTT run. For each of the comets with multi-run data the individual results are presented here; the conclusions 
based on the combined data, including phase curves, are given in the text. 

^ Faint coma present. Note that these results do not include any correction for the presence of faint near-nucleus coma. 

e Results from 2006 INT run. 

/ Results from 2006 FTN observations. 
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Figure 15. Image showing 36P, made up of 17x200 second expo- 
sures taken on the 6th March 2005. There is some possible activity 
in the image, which appears to be confirmed by the active profile, 
where a clear difference between the comet and background PSF 
is visible. 



Rh > 4 AU is surprising for a comet that has never been 
seen to show any large amount of activity. 

It is clear that the total flux is still dominated by the 
flux from the nucleus, as obvious variations can be seen in 
the brightness of the comet in fig. 1181 which shows the vary- 
ing brightness in the original data folded onto a period of 
3.5 hours. This period corresponds to twice the fitted pe- 
riod of 0.074 days given by the strongest minimum in the 
periodogram (fig. I16[) ; a minimum of y 2 jv — 0.55, which is 
actually below the expected value, but within 2a of the ex- 
pected x 2 l v — !• The null hypothesis of constant brightness 
is rejected at a 6<r level, implying that the variations are 
real. 

The unknown contribution to the flux from the coma 
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Figure 16. Periodogram for 36P, using the total observed bright- 
ness of the comet. 

causes uncertainty in both the average magnitude of the 
comet and the amplitude of the light-curve variation. As- 
suming a steady state coma, its contribution is modelled 
as a constant flux level, which is then subtracted from the 
data to leave the variation due to the nucleus. A constant 
coma contribution acts to raise the overall brightness, and 
to reduce the amplitude of the observed variation, as it will 
contribute relatively more when the nucleus cross section 
is at a minimum than at maximum. Taking the measured 
coma contribution from the surface brightness profile (40% 
of the average total flux), and removing this from the data, 
gives approximate 'nuclear' magnitudes. Figures [T7l and [181 
show the periodogram and folded light-curve for these data. 
The rotation period found does not change; the presence or 
lack of coma at this level does not affect the measurement 
of this value. It can be seen that the amplitude of the light- 
curve is greatly increased from Am — 0.4 mag to Am = 
0.65 mag, and that this increase in the size of variation rel- 



16 C. Snodgrass et al. 




o I ' ' ' ' 1 ' ' ' ' 1 ' ' ' ' 1 

0.1 0.2 0.3 

Period (days) 

Figure 17. Periodogram for 36P, from data with a constant flux 
level due to the coma removed. 
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Figure 18. 36P photometric data folded onto a 3.5 hour rotation 
period. The lower data set, with the dashed line drawn through 
it, are the data with a constant coma flux equivalent to m c = 
22.4 removed from them. 



ative to the size of the error bars on individual data points 
increases the X 2 /v- There are no longer any minima within 
X 2 /v = 1 ± y/2/v; the strongest minimum in the coma sub- 
tracted periodogram has \ 7 ' l v ~ 1-8, corresponding to 2.3<r. 
The mean magnitude is increased by removing a coma com- 
ponent; the original data has tur — 21.37, removing 40% of 
this flux gives ma — 21.94. These values are all calculated 
taking the value of 40% to be fixed; in reality the 23% un- 
certainty on the flux contribution leads to large error bars 
on any coma subtracted photometry. 

The sizes and shapes implied by the two light-curves 
are different. The original data give a radius for the equiv- 
alent sphere of 2.5 km, and an axial ratio of a/b > 1.4, 
and therefore nuclear dimensions of a x b = 3.3 x 2.3 km. 
The coma subtracted model gives a smaller, more elongated 
body, with tn = 1.9 km, a/b > 1.8 and thus a x b — 3.1 x 1.6 
km. The first solution matches the va lue of tn = 2.5±0.2 km 
found bv lLowrv fc Weissm an (2003), who obtained a snap- 
shot of 36P when it was inbound at 4.4 AU in May 2001. 
The second solution is also consistent with this considering 
the uncertainties on the coma fraction. 



Figure [T8l shows that the faint coma detected does not 
prevent measurement of a rotation period, and that this pe- 
riod is independent of any corrections to the photometry 
applied to correct for coma. It is the determined period 
which is of greatest interest for this object; at 3.5 hours, 
these data show 36P to be the fastest rotating comet known. 
This implies that it has a surprisingly dense nucleus, under 
the assumption that it is essentially strengthless, or that 
this nucleus has strength not seen in others. Either of these 
present a challenge to the general picture of nuclei that is 
beginning to be built up; that they are weak and very low 
density bodies. The required density for a strengthless body 
rotating at a 3.5 hr period depends on its elongation, and 
therefore for 36P depends on whether or not a coma model 
is subtracted. Equation gives minimum bulk densities for 
36P of 1.2 g cm" 3 for the original data (with a/b > 1.4) 
and 1.6 g cm" 3 for the coma subtracted version. Clearly the 
density is greater than water ice, for any level of coma, and 
over twice that required for any other nucleus. 

Colours were also measured for 36P, and were found 
to be consistent across both nights with average values of 
(V-R) = 0.46±0.05 and (R-I) = 0.52±0.04. These colours 
are entirely typical of JFC nuclei and fall near the centre of 
the observed distribution, although it must be remembered 
that some fraction of this colour will be due to dust grains 
in the detected coma. 

5.2 INT data: March 2006 

With this potentially interesting result in mind, 36P was 
the primary target of the 2006 INT run. After the first two 
nights were lost to poor weather, 36P was intensively mon- 
itored on the third to ensure unambiguous detection of any 
3.5 hour period. A total of 63 r'-band images were taken 
over ~ 5.5 hours, meaning that any short period would 
give almost complete phase coverage. As this night was non- 
photometric, and to search for longer periods, a considerable 
amount of data were also gathered on 36P on the fourth 
night, with 55 r'-band frames taken over almost 7 hours. A 
combination of these data and a surface brightness profile 
of the resultant image (fig. I19p show that the comet was 
effectively a bare nucleus; a year later and 0.7 AU further 
from the Sun the faint activity seen in the NTT data had 
ceased. The formal coma limit measured within p — 5" is 
m c > 23.31, or < 20% of the flux, however this is probably 
an over-estimate due to the slight rise in the outer part of 
the profile due to residual sky noise (note that the profile 
is entirely consistent with the stellar PSF within the error 
bars) . 

Although more distant during the INT run than in the 
NTT run, 36P was closer to opposition, and consequently of 
similar apparent magnitude: ran. = 21.570±0.008. Assuming 
the standard linear phase law with j3 = 0.035 mag. deg -1 
gives ma (1,1,0) = 15.247 for the bare nucleus; compari- 
son with the non-coma subtracted value from the NTT run, 
771^(1,1,0) = 15.272, allows calculation of the fraction of 
flux from the nucleus during that run by: 

^Q0.4(m R ~m N ) _ ^Q0.4(m n (l,l,0) NTT -m R (l,l,0) INT ) _ ^Q2% (6) 

Clearly this cannot be true: the bare nucleus should be 
fainter than the active one, and the fraction of flux due to 
the nucleus should be ~ 60% if the measurement of 40% 
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Figure 19. Image showing 36P, made up of 55x75 second ex- 
posures taken on the 2nd March 2006 with the INT, and the 
corresponding profile. The comet is clearly inactive. 
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Figure 20. Periodogram for 36P from INT data. 
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Figure 21. 36P differential photometric data folded onto an 40.6 
hour rotation period. 
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Figure 22. 36P differential photometric data taken on the 1st 
March 2006. 
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Figure 23. 36P differential photometric data taken on the 2nd 
March 2006. 



coma contamination is correct. The discrepancy is probably 
due to an incorrect phase law. A discussion on the phase 
function for 36P is given in the next section. 

The full range of the INT data is Am = 0.7 ± 0.1 mag, 
larger than that measured in the NTT data but similar to 
the range observed in the coma-corrected data. It implies 
a/b > 1.9. The periodogram for the INT data is shown in 
fig. 1201 which shows that the fit to the fluctuations in the 
brightness is dominated by a long period variation. This pe- 
riodicity search is based on differential magnitudes, since the 
large field of view of the WFC meant that the same stars 
could be used for comparison on both nights, meaning that 
no uncertainty is added due to calibration onto a standard 
scale. The data folded onto the best period, at P rot = 40.6 
hours, is shown in fig.[3T] Such a long period could not have 
been detected in the shorter and sparser data set obtained 
with the NTT, but there remains the question of where the 
3.5 hour period came from, and whether it has any basis 
in reality. It is clear from fig. [16] that there are many sta- 
tistically acceptable periods in the NTT data; due to the 
few data points many periods can be found which have suf- 
ficiently low X*/v. Tne NTT data can be folded onto a 40.6 
hour period to give a visually acceptable light-curve, how- 
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Figure 24. Image showing 36P, made up of 84x120 second ex- 
posures taken on the 26th February 2007 with the NTT, and the 
corresponding profile. The comet is clearly inactive. 
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time (MJD 54150+ ) 

Figure 25. 36P differential photometric data taken on the 26th 
February 2007. The slight rise at i sj 8.29 is due to proximity 
to a faint background star. Overall there is a drop of about 0.2 
mag. over the ~ 4 hours. 



ever this could be said of a large number of periods given 
the sparse data. Of more interest is whether or not any short 
period variation exists in the INT data. The unfolded data 
from each of the INT nights are shown in fig. [22] and fig. 1231 
which show hints of shorter period variation. It is possible 
that the 3.5 hour period and these short variations are part 
of a non-symmetric light-curve with a longer period. The 
relatively noisy INT data does not allow a definite conclu- 
sion that the 40.6 hour period is correct, nor entirely rule 
out a short period. 

A large number of colour frames were taken within the 
time-series on 36P, giving good phase coverage for the fit- 
ted (long) period. Using the calibrated magnitudes of stars 
from the fourth (photometric) night gave the colours from 
differential magnitudes involving the same stars on the non- 
photometric night. All individual colours are consistent with 
means of (V - R) = 0.48 ± 0.03 and (7? - J) = 0.62 ± 0.02, 
within the error bars on individual points. 



5.3 NTT data: February 2007 

One of the main conclusions of this paper must unfortu- 
nately be that it is simply not practical to reach the re- 
quired S/N for accurate nucleus light-curves with a 2.5m 
telescope. We therefore decided to return to the 3.6m NTT 
to make a final search for short period variation. We were 
awarded 4 hours of Director's Discretionary Time (DDT) 
and monitored 36P continuously over 3.7 hours, and also 
returned to it for 15 minutes on the following night. Con- 
ditions were excellent and the comet was easily detected 
with m_R = 22.774 ± 0.013, with the comet near aphelion at 
Rh = 5.2 AU. The combined image and profile appear inac- 
tive (fig. [24}, although the profile is noisy in the outer part 
at E ~ 30 mag arcsec -2 , and the formal limit is m c > 24±2, 
or < 22 ± 50% of the flux. 

There was no short period variation found in the dif- 
ferential light-curve (fig. [25}. The short sequence of data 
on the second night was taken to search for longer period 
variations, but the light-curve remained fairly flat over both 
nights. Statistically, the null hypothesis of a constant mag- 
nitude with no variation is entirely acceptable: x 2 /v = 1.38, 
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Figure 26. 36P mean magnitudes reduced to rg = A = 1 
AU and plotted as a function of phase angle a. The labels 
give the telescope and year of the time-series. The two la- 
belle d 'NTT ss' are the 2007 sn ap-shots, the point 'LW03 ss' is 
fromlLo wrv &; Weissmanl l l2003h and the point 'E+07 ss' is from 
lEpifani et al .1 (|2007t >. with the arrow indication that the comet 
was faintly active when observed by these authors and therefore 
this point is a lower limit. Each error bar is the combination of 
the photometric uncertainty on the point and the uncertainty due 
to only observing a section of the rotational light-curve. The error 
bar on the NTT 2005 point includes the uncertainty on the coma 
flux, which the point has been corrected for. 



within the 1 ± 0.15 range on a la variance, however there is 
a linear trend across the 3.7 hours of the first night, with a 
drop of Am w 0.2 mag. We suspect that this is a short seg- 
ment of a much longer period light-curve. 36P rotates with 
a long period; although we still cannot be sure of what this 
period is, it appears to be > 24 hours. 

The (V — R) colour of the nucleus was measured to 
enable calibration of the photometry including the colour 
term: From a single U-band snap-shot in the middle of the 
7?-band time-series we obtain (V — R) = 0.40 ± 0.17. A 
weighted average of all colour measurements gives (V — R) = 
0.47 ± 0.02 and (R - /) = 0.60 ± 0.02. 
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The three independent light-curve measurements do al- 
low the calculation of an approximate phase function. We 
take the three measurements (2005 NTT at a = 11.1°, 2006 
INT at a = 1 .0°, 2007 NTT at a = 6.7° ), together with the 



snap-shot bv lLowrv fe Weissman ( 20031) taken at a = 7.3°, 
the deep image from lEpifani et af 1 20071 ) (when the comet 
was weakly active at Rh — 3.9 AU and a = 9.3°) and two 
further snap-shots taken by the authors with the NTT dur- 
ing 2007 (at a = 10.7 and 11.0°, both at R h > 5 AU and ap- 
parently inactive). Figure [26] shows this data as mn(l, 1, a) 
against a. The error bars for each point cover the full range 
of uncertainty due to both the photometric uncertainty (very 
low for the light-curves as it decreases as %/iV) and also the 
uncertainty due to the unknown rotational phase of the ob- 
servations. The maximum Am observed is 0.7 mag., so the 
snap-shots can be up to 0.35 mag. from the mean magnitude 
and the DDT partial light-curve covering a 0.2 mag. range 
can be up to 0.25 mag. from the mean. The 2005 NTT data 
point is given with a 40±23% coma contribution subtracted, 
and the error bars on this point include this considerable un- 
certainty. Note that we believe the INT data point to be well 
determined, as reflected by the small uncertainties on it, as 
the large amount of data over two nights covered Am = 0.7 
and thus gives a good average magnitude. 

The best fit linear phase function has = 0.060 ± 0.019 
mag. deg -1 . This is quite steep, matching the previous 
largest reported value fo r JFCs of /3 = 0.06 for 2P/Encke 
(Fernande z fc et al.ll2000l ). We fit only a linear phase func- 
tion as the data do not justify higher order fits, s uch as 
the standard IAU (H, G) system (jBowell et al.lll989h . Also, 
while the INT data point is at low phase angle, we do 
not expect the sort of opposition surge seen in asteroids; 
iRabinowitz et all l|2007h show that for outer Solar System 
bodies (TNOs and Jupiter Trojans) there is no opposition 
surge. 

The size of the nucleus is determined from the absolute 
magnitudes from this phase law; extrapolating the linear 
phase laws back to a — 0° gives rna{l, 1, 0) = 15.223±0.008. 
Assuming a 4% albedo this gives a radius of 2.58 ± 0.01 
km. When taken with a/b = 1.9 this implies dimensions of 
a x b = 4.1 x 2.1 km. 



6 DISCUSSION OF COMET PROPERTIES 
AND COMPARISON WITH KBOS 

Here we discuss the effect of this new data on the general 
trends for JFCs id entified in paper II. In that paper, and in 
lLowrv et al] (§003) we suggested that the cut off in P ro t-a/b 
space, and corresponding cut off in densities, corresponded 
to the true density of a typical JFC nucleus, at ~ 0.6 g 
cm" 3 . If 36P was rotating in only 3.5 hours, it challenged 
this theory, however we now believe that it is one of the 
slower rotators and does not help to define the cut off. 94P 
also appears to be a slow rotator. The possible periods for 
121P all place it well within the previously observed range, 
while the lack of period determinations for the other comets 
means that our conclusions on the rotational properties of 
JFC nuclei are unchanged from paper II. 

The new results presented here do allow us to update 
the colour-colour plots presented in the earlier papers. The 
trend of increasing (R — I) with increasing (V— R) described 
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Figure 27. (R — I) against (V — R) for all JFCs with known 
colours. Filled circles — this work; open circles - comets with 
colours previously determined using the same multi-filter pho- 
tometry method used here (papers I & II and references therein). 
The position of the Sun on these axes is marked using the symbol 
0. 



in these papers is not as visually obvious with the inclusion 
of the INT data (fig. I27[) . however this is largely due to the 
weakly active comets with large error bars on the colours. 
A weighted best fit to all the data still gives a straight line 
consistent with our earlier result at a 2a level: 

(t, r, , I7 m f m = 0.74 ±0.12 

(R-IU c = m(V-R) mc + ci c = 012±001 (7) 

While comets 45P/Honda-Mrkos-Pajdusakova, 
P/2004 H2 and 40P are not very close to this line 
(the first two well below the line in the lower centre of 
the figure, 40P above the trend in the upper left) they 
all showed some activity during the observations. In the 
case of both 40P and P/2004 H2 the coma was very weak, 
while 45P was observed by iLamv et al.1 (l999i) using the 
HST when the comet was active and at large phase angle 
(a « 90°), but with a coma subtraction method. In all of 
these cases, the activity is not sufficient to reject the colour 
data as weak activity at large Rh should be dust dominated 
and therefore of similar optical properties to the nucleus, 
but also means that their difference from the rest of the 
data should not be over-interpreted. The mean colours of 
JFC nuclei are (V - R) nuc = 0.45 ± 0.11 (N = 31) and 
(R - I) nuc = 0.46 ± 0.10 (N = 19). 

The updated fit remains consistent with the equivalent 
fit to KBO data (fig. gSJ, while it is clear that the KBOs 
have redder colours. To further investigate the idea that the 
more neutral JFCs are 'resurfaced' KBOs, it is desirable to 
use a one dimensional parameter to describe the relative 
redness of these objects. To do this the (V — R) and (R — I) 
data shown in fig. [2H] were translated into distance along the 
best fit line x, which has arbitrary units linearly related to 
the normal colour indices, with increasing x corresponding 
to redder colours. The weighted average best fit between the 
JFC and KBO best fit lines quoted above was used (with 
m = 0.90 ± 0.05 and c = 0.09 ± 0.01) giving a translation of 



x = 0.45(1/ -R) + 0.67 (R -I) - 0.06 



(8) 
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Figure 28. (R - /) against (V - R) for all JFCs, KBOs and 
Centaurs with known colours. Open circles - JFCs; filled circles - 
KBOs; triangles — Centaurs. The position of the Sun on these axes 
is marked using the symbol 0. The dashed line shows our best fit 
to the comet data, the dotted line i s the fi t to the KBOs. KBO an d 
Centaur data from ljewitt fc Luul feOOlft ; [Peixinho et al.l i2004) . 




more neutrally coloured objects remaining. A least squares 
fit gave f(x) — 2.8 — 3.5a; (again, in non-physical units) and 
gave a good reproduction of the observed JFC distribution. 
This de-reddening function means that 100% of the KBO 
surfaces remain at x — 0.5 and they are entirely depleted 
beyond x = 0.8. Clearly this model is very approximate in 
its treatment of distributions, but it serves to demonstrate 
a possible transformation route between the KBO and JFC 
distributions. A more thorough approach would be to treat 
a large number of test particles, with a starting distribution 
like that observed in KBO colours, and then test the effect of 
different theoretical 'de-reddening' models on the resultant 
JFC colours by computing the change in colour for individ- 
ual comets. While an interesting problem, such calculations 
are beyond the scope of this work. Further colour data, es- 
pecially on the JFCs, are required before a more detailed 
analysis of the distribution is justified. 

Although a linear decrease in the number of surviving 
KBO surfaces with colour can be invoked to explain the ob- 
served JFC colours, a physical interpretation of this effect 
is required. A plausible solution is to assume that the ob- 
served variation in KBO colour is associated with the ages 
of their surfaces, and that space weathering reddens the sur- 
faces in an approximately linear fashion with time. Under 
this paradigm the objects in very red tail of the distribution 
are very old, while the bluer ones have under gone some 
'resetting' of their surfaces (e. g. due to a collision) more 
recently. When a KBO is perturbed into the inner Solar 
System, the onset of cometary activity acts to resurface the 
nucleus. For younger KBOs, whose surfaces are not very dif- 
ferent to their sub-surfaces as they have not been weathered 
for long, this has little effect, but for the very ancient and red 
surfaces, the difference is important. One interesting aspect 
of this interpretation is that it suggests that nuclei surfaces, 
being younger, are a closer match to their interiors, and thus 
may provide a better view of 'primordial' material than the 
ancient weathered surfaces of KBOs. 



Figure 29. The (R — /) against (V — R) data were reduced onto 
distance x along the best fit (see text). This plot shows the his- 
tograms of N(x) against x for JFCs and KBOs. The KBOs (clear) 
are reasonably well matched by a normal distribution (dashed 
line). The JFCs (hatched) appear to have a truncated distribu- 
tion; the model JFC population (solid line) was created by con- 
volving the KBO model with a 'de-reddening' term that is oc —x. 

Histograms of the number of objects N(x) per bin (dx = 
0.05) are shown in fig. [29] It can be seen that the JFC dis- 
tribution does indeed appear to be consistent with the KBO 
distribution, but lacking in 'ultra-red' matter. The two pop- 
ulations are a good match at the bluer end. The KBO colours 
appear to be approximately normally distributed (the fitted 
distribution has mean 37kbo = 0.75 and standard deviation 
ctkbo = 0.14), as do the JFCs, despite the relatively poor 
statistics (xjfc = 0.61, ctjfc = 0.10). The model distri- 
bution plotted for the JFCs (solid line) however is not the 
best fit normal distribution, but one generated by convolv- 
ing the fitted KBO distribution with a 'de-reddening' term: 
N(x)jfc = N (x) kbo <Sif{x), where the fraction of objects re- 
maining f(x) oc —x. Thus the 'de-reddening' depends on the 
colour, with the most red KBOs being most depleted, with 



7 SUMMARY 

We report here observations on JFCs from two runs with 
the Isaac Newton Telescope, with additional data from the 
Faulkes Telescope North. A large number of snap-shots 
were taken, and in addition time-series data were taken 
on 5 comets. We also describe our observations of comet 
36P/ Whipple over the last 3 years, at both the INT and 
ESO's NTT. Our results are as follows: 

(i) For the majority of the comets we have snap-shot data. 
Of these, five of the comets were not detected, three were 
visibly active, and the remainder (9 comets and one Damo- 
cloid) were apparently inactive at large distance. Most of 
these results were close to the detection limits of the 2.5m 
INT; we must conclude that good S/N nucleus observations 
really require a 4m class telescope. The magnitudes, im- 
plied sizes, colours and activity levels for each comet are 
summarised in tabled 

(ii) For the comets which were observed over a time- 
series, constraints were also placed on the rotation period 
and the elongation of the nucleus, from the range in magni- 
tudes observed. Faint activity around 40P and P/2004 H2 
prevented a measurement of these parameters for these 
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comets, but for 47P, 94P and 121P and we measure min- 
imum axial ratios of a/6 > 1.4,3 and 1.1 respectively. This 
implies that 94P has a very elongated nucleus. Periodogram 
analysis found best fit rotation periods of ~ 33 and 10 hours 
for 94P and 121P, but both of these are very uncertain and 
should be regarded as estimates. 

(iii) For 36P we present data taken at 5 epochs: 3 (par- 
tial) light-curves and 2 snap-shots. The first data set gave 
a preliminary result implying fast rotation for this nucleus, 
however larger and better sampled later data sets revealed 
that the rotation period is long (> 24 hours). It is still not 
well determined; the largest data set has a best fit period 
of ~ 40 hours, but this is not a unique solution. This data 
set also shows a large amplitude variation implying an axial 
ratio a/b > 1.9. 

(iv) The data on 36P were obtained at a range of phase 
angles 1 < a < 11°. Including the uncertainty on the snap- 
shot/short time series data due to the incomplete phase cov- 
erage, all the points are consistent with a linear phase func- 
tion with p — 0.060 ± 0.019 mag. deg~ . We also measure 
an approximate phase function for 47P, with a very steep 
(i = 0.083 ± 0.006 mag. deg" 1 . 

(v) U-RJ-band photometry was obtained on seven of the 
comets, including 36P, and the distribution of their colours 
was investigated. JFC nuclei and KBOs are found to have a 
similar trend of increasing (R — I) with increasing (V — R), 
with the JFCs being less red in general. We define a single 
parameter x to measure the redness of the objects across the 
VRI bands. We show that the distribution of JFC colours 
provides a good match to the blue end of the distribution 
of KBO colours, and that it is possible to recreate the JFC 
distribution by applying a 'de-reddening' function that goes 
as f(x) oc —x to the KBO distribution. One interpretation 
of this model is that it relates to the age of the comet. 
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